Draft version February 26, 2008 

Preprint typeset using I^T^X style cmulateapj v. 10/09/06 



oo 
O 
O 
(N 

X> 

PL, 
VO 

Or 

6 



BIMODALITY IN DAMPED Lya SYSTEMS 
Arthur M. Wolfe 1,2 , Jason X. Prochaska 1 ' 3 , Regina A. Jorgenson 1,2 & Marc Rafelski 1,2 

Draft version February 26, 2008 

ABSTRACT 

We report evidence for a bimodality in damped Lya systems (DLAs). Using [C II] 158 /im cooling 
rates, £ c , we find a distribution with peaks at £ C =10~ 27A and 10 -26 ' 6 ergs s _1 H _1 separated by a 
trough at £" xt « 10~ 270 ergs s _1 H _1 . We divide the sample into 'low cool' DLAs with £ c < £^ Tlt 
and 'high cool' DLAs with £ c > ££ rlt and find the Kolmogorv-Smirnov probabilities that velocity 
width, metallicity, dust-to-gas ratio, and Si II equivalent width in the two subsamples are drawn from 
the same parent population are small. All these quantities are significantly larger in the 'high cool' 
population, while the H I column densities arc indistinguishable in the two populations. We find that 
heating by X-ray and FUV background radiation is insufficient to balance the cooling rates of either 
population. Rather, the DLA gas is heated by local radiation fields. The rare appearance of faint, 
extended objects in the Hubble Ultra Deep Field rules out in situ star formation as the dominant 
star-formation mode for the 'high cool' population, but is compatible with in situ star formation as 
the dominant mode for the 'low cool' population. Star formation in the 'high cool' DLAs likely arises 
in Lyman Break galaxies. We investigate whether these properties of DLAs are analogous to the 
bimodal properties of nearby galaxies. Using Si II equivalent width as a mass indicator, we construct 
bivariate distributions of metallicity, £ c , and arcal SFR versus the mass indicators. Tentative evidence 
is found for correlations and parallel sequences, which suggest similarities between DLAs and nearby 
galaxies. We suggest that the transition-mass model provides a plausible scenario for the bimodality 
we have found. As a result, the bimodality in current galaxies may have originated in DLAs. 
Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 
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1. INTRODUCTION 

Recent studies of ~ 10 5 low-redshift galaxies reveal a 
striking bimodality in their properties. The galaxy pop- 
ulation is divided into two mass sequences; a 'blue' se- 
quence with stellar masses M* < 10 10 5 M Q and a 'red' 
sequence in which M* > 10 10 5 M Q (Kauffmann et al. 
2003; Baldly et al. 2004). The 'blue' sequence is com- 
prised of late-type galaxies undergoing active star for- 
mation, while the 'red' sequence consists of early-type 
galaxies with little, if any, star formation. Dekel & Birn- 
boim (2006) suggest that the mass sequences reflect the 
different paths for the build up of stellar populations 
in galaxies. The critical parameter in their scenario is 
the dark-matter halo 'shock mass' , M s i 1O ck=10 11 ' 5 M , 
which corresponds to M* = 10 10 5 Mq. At high redshifts, 
gas accreted onto halos with masses -Mdm < Mshock is 
not heated to the virial temperature of the halo, but 
rather produces star-forming disks through the infall of 
cold streams. The continuation of this process to the 
present results in the formation of the 'blue' sequence. 
By contrast gas accreted onto halos with masses Mdm 
> M s hock is shock heated to the virial temperature, and 
accretes onto the halos in a hot cooling flow. Cold fila- 
ments penetrate the hot gas and their inflow results in 
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star-forming bulges in massive galaxies (Dekel & Birn- 
boim 2006; see also Bell et al. 2004; Faber et al. 2007). 
In this case the low density of the hot gas and feedback 
processes suppress star formation at z < 2. As a result 
these halos evolve onto the 'red' sequence. 

Our purpose here is to present evidence for an analo- 
gous bimodality in the properties of damped Lya sys- 
tems (DLAs). Because these neutral-gas layers are the 
likely progenitors of modern galaxies (see Wolfe, Gawiser, 
& Prochaska 2005 [hereafter referred to as WGP05] for 
a review), the DLA bimodality may be related to the 
galaxy phenomenon discussed above. In fact we shall ar- 
gue that the bimodality in modern galaxies originates in 
DLAs. We use a technique which, for the first time, mea- 
sures bimodality in absorption- line gas 4 . Rather than 
obtain signatures based on starlight emitted by galaxies 
(e.g. Kauffmann et al. 2003), we rely on an absorption- 
line diagnostic, which we argue is a signature of star 
formation. Specifically, we measure the C II* A1335.7 
absorption line, which, if thermal balance is assumed, is 
an indicator of the rate at which neutral gas is heated. 
The crucial parameter is the [C II] 158 fim cooling rate 
per H atom, £ c , which divides the DLA sample in such a 
way that objects with £ c > £™' differ fundamentally from 
those with 4 < £ c / lt , where £ c c rit -10" 27 ergs s" 1 H _1 . 
We shall argue that £ c is a tracer of star formation rates 
(see Wolfe, Prochaska, & Gawiser 2003 [hereafter referred 
to as WPG03]) and that the bimodality in £ c is physically 
related to a transition in star-formation modes, from in 
situ star formation in DLAs with £ c < £" lt to star for- 

4 Prochaska et al. (2002) report evidence for bimodality in the 
ratio of the nitrogen to alpha element abundance, but its statistical 
significance is tentative owing to the small size of the sample. 
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mation in compact 'bulge' regions sequestered away from 
the DLA gas in objects in which £ c > £" lt : the higher star 
formation rates (SFRs) predicted for the 'bulge' models 
account for the higher heating rates of the surrounding 
DLA gas (see § 7) . We shall argue that the bimodality in 
star formation modes is caused by a transition in galaxy- 
formation modes, which ultimately is due to a transition 
in mass (Dekel & Birnboim (2006). 

The paper is organized as follows. In § 2 we examine 
the £ c distribution and discuss the results of tests to as- 
sess whether the distribution is bimodal. In § 3 we divide 
our DLA sample into two sub-samples: 'low cool' DLAs 
with 4 < £ c c rit and < high cool > DLAs with £ c > ^crit_ For 

each sub-sample we compile distributions of physical pa- 
rameters such as absorption-line velocity width, metallic- 
ity, dust-to-gas ratio, Si II A 1526 equivalent width, and 
H I column density Am- We then determine the prob- 
ability that the pair of distributions corresponding to a 
given parameter is drawn from the same parent popu- 
lation, and describe the results of these tests. In § 4 
we discuss the physical significance of £™' and conclude 
that it signifies a dividing line between two modes of 
star formation. In § 5 we give a brief summary of the 
results and conclude that the crucial parameter distin- 
guishing the two populations is dark-matter mass. In § 
6 we draw analogies between bimodality in DLAs and in 
modern galaxies and show how the two are related. We 
also place these results in the context of modern theories 
of galaxy formation. Conclusions are given in § 7. 

Throughout this paper we adopt a cosmology with 
(fi M ,n A ,fc)=(0.3,0.7.0.7). 
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Fig. 1. — £ c versus TVhi for the 76 DLAs in our sample. Red data 
points are positive detections, green are upper limits, and blue are 
lower limits. The blue star is the average value for the Galaxy disk. 



2. BIMODALITY OF THE l c DISTRIBUTION 

In this section we describe evidence for bimodality in 
the distribution of 158 /mi cooling rates, l c . The values 
of £ c for the 76 DLAs in our sample are listed in Table 
1 (column 7). As in previous papers (e.g. WPG03) we 
define £ r as follows: 
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Fig. 2. — Histogram depicting distribution of the 37 positive 
detections of l c reported in Table 1. 



where N(C II*) is the column density of the excited 2 P^/2 
state in the 2s 2 2p term of C + , A u i is the Einstein coeffi- 
cient for spontaneous photon decay to the ground 2 P\/2 
state (A u i=2Ax 10 -6 s _1 ), and hv u i is the energy of this 
transition (hv ul /k=92 K). All of the N(C n*) values in 
Table 1 (column 6) were deduced from velocity profiles 
of the C II* A 1335.7 transition arising from the 2 -p3/2 
state. We obtained 47 profiles with the High Resolution 
Echelle Spectrometer (HIRES; Vogt et al. 1994) on the 
Keck I 10 m telescope. The HIRES spectra were ac- 
quired using either a 0.8 " or 1.1 " wide decker, resulting 
in velocity resolution with FWHM = 6 and 8 km s _1 
respectively. Sixteen of the profiles were obtained with 
the Echellette Spectrograph Imager (ESI; Sheinis et al. 
2002) on Keck II. The ESI spectra were obtained with 
a 0.5 " or 0.75 " slit (FWHM w 33 and 44 km s~ x re- 
spectively). The remaining 13 profiles were obtained by 
other observers with the UVES spectrograph on the VLT 
8 m telescope. The Am values were derived by fitting 
Voigt profiles to damped Lya lines detected in the SDSS 
survey (Prochaska et al. 2005) and with the Keck tele- 
scopes, the VLT, and several 4 m class telescopes (Peroux 
et al. 2003; Storrie-Lombardi & Wolfe 2000; Wolfe et al. 
1995). The sample in Table 1 comprises 32 upper lim- 
its (95 % confidence level), 37 positive detections, and 7 
lower limits on £ c . 

We first noticed evidence for bimodality in previous 
plots of £ c versus Ahi (e.g. WGP05). An updated ver- 
sion of these plots is shown in Fig. |TJ The red data points 
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Q0405-443 


00:00:00.00 


00:00:00.00 


2.5950 


20.901^ 


13.66 ± 0.21 


-26 


.76 + 0.23 


79 


-0.96 + 0.10 


-1.33 + 0.10 




24 


Q2359-02 


00:01:50.00 


-01:59:40.34 


2.1539 


20.301°;!° 


< 14.48 


< 


-25.34 


78 


-1.58 + 0.01 


-1.88 + 0.03 


0.28 


6,13 


Q2359-02 


00:01:50.00 


-01:59:40.34 


2.0951 


20.70!°;!° 


13.70 ± 0.06 


-26 


51 + 0.12 


142 


-0.77 + 0.02 


-1.65 + 0.03 


0.83 


6,13 


J001328. 21+135827 


00:13:28.21 


+13:58:27.9 


3.2811 


21.55+°;!* 


13.67 + 0.04 


-27. 


.40 + 0.16 


10 


-2.10 + 0.16 


-2.72 + 0.02 


0.28 


32 


BR0019-15 


00:22:08.01 


-15:05:38.78 


3.4389 


20.92+°;!°, 


13.84 + 0.02 


-26 


.60 + 0.10 


118 


-1.06 + 0.05 


-1.58 + 0.04 


0.76 


6,13 


PH957 


01:03:11.38 


+13:16:16.7 


2.3090 


91 ,-+0.08 

-0.08 


13.59 ± 0.05 


-27 


.30 + 0.09 


56 


-1.46 + 0.01 


-1.90 + 0.04 


0.38 


1,6,13 


SDSS0127-00 


01:27:00.69 


-00:45:59 


3.7274 


21.15+°;!° 


13.20 + 0.06 


-27 


46 + 0.12 


40 


-2.40 + 0.10 


-2.90 + 0.02 


0.28 


22 


PSS0133+0400 


01:33:40.4 


+04:00:59 


3.7736 


20.55+°;!° 


14.02 ± 0.01 


-26. 


.05 + 0.10 


123 


-0.75 + 0.10 


-1.07 + 0.08 


0.92 


22 


PSS0133+0400 


01:33:40.4 


+04:00:59 


3.6919 


20.701°;!° 


12.95 ± 0.03 


-27 


.27 + 0.11 


39 


-2.34 + 0.15 


-2.74 + 0.05 


0.31 


22 


J013901. 40-082443 


01:39:01.40 


-08:24:43.9 


2.6773 


20.70+°!g 


13.81 ± 0.03 


-26. 


.40 + 0.15 


110 


-1.27 + 0.19 


-1.62 + 0.02 


0.67 


32 


Q0149+33 


01:52:34.472 


+33:50:33.23 


2.1408 


20.50+°,;!°, 


< 12.78 


< 


-27.24 


40 


-1.49 + 0.05 


-1.77 + 0.02 


0.22 


6,13 


Q0201+11 


02:03:46.53 


11:34:40.4 


3.3869 


21 26+ ' 10 


14.12 ± 0.10 


-26. 


.66 + 0.14 


67 


-1.25 + 0.11 


-1.41 + 0.05 


0.64 


15 


PSS0209+0517 


02:09:44.52 


+05:17:17.3 


3.8636 


20.55_ 1() 


< 12.51 


< 


-27.55 


47 


-2.60 + 0.10 


-2.96 + 0.09 


0.06 


22 


SDSS0225+0054 


02:25:54.85 


+00:54:51 


2.7137 


21.001°;!* 


> 13.37 


> 


-27.15 


60 


-0.91 + 0.14 


-1.31 + 0.04 


1.08 


29 


J023408.97-075107 


02:34:08.97 


-07:51:07.6 


2.3180 


20.951°;!* 


< 13.41 


< 


-27.05 


10 


-2.74 + 0.14 


-3.14 + 0.04 


0.09 


32 


J0255+00 


02:55:18.62 


+00:48:47.94 


3.9146 


21 30+ - 05 


13.44 + 0.04 


-27. 


.37 + 0.06 


38 


-1.78 + 0.01 


-2.05 + 0.09 


— 


13 


J0307-4945 


03:07:22.85 


-49:45:47.6 


4.4679 


20 67+ - 09 


< 13.59 


< 


-26.60 


192 


-1.55 + 0.08 


-1.96 + 0.21 


— 


11 


Q0336-01 


03:39:00.99 


-01:33:18.07 


3.0621 


21 20+ ' 10 


14.01 ± 0.01 


-26. 


.71 + 0.10 


102 


-1.54 + 0.01 


-1.81 + 0.02 


— 


13 


Q0347-38 


03:49:43.54 


-38:10:04.91 


3.0247 


20 63+° °° 


13.47 + 0.03 


-26 


68 + 0.03 


84 


-1.16 + 0.03 


-1.62 + 0.01 


0.45 


6,13,20 


Q0458-02 


05:01:12.81 


-01:59:14.25 


2.0396 


2i.65l°;° 9 9 


> 14.80 


> 


-26.37 


82 


-1.19 + 0.02 


-1.76 + 0.05 


0.67 


6,13 


HS0741+4741 


07:45:21.75 


+47:34:35.56 


3.0174 


20.48+°,;!°, 


< 12.55 


< 


-27.44 


42 


-1.68 + 0.00 


-1.93 + 0.00 


0.21 


13 


FJ0747+2739 


07:47:11.19 


+27:39:03.6 


3.9000 


20.501°;!° 


13.35 ± 0.07 


-26. 


.66 + 0.12 


150 


-2.01 + 0.01 


-2.45 + 0.03 


0.22 


22 


FJ0812+32 


08:12:40.8 


+32:08:08 


2.6263 


21.351°;!° 


14.30 + 0.01 


-26 


56 + 0.10 


70 


-0.93 + 0.05 


-1.76 + 0.01 


0.60 


22,31, 


J081435. 18+502946 


08:14:35.18 


+50:29:46.5 


3.7082 


21.351°;!* 


< 13.18 


< 


-27.69 


30 


-3.00 + 0.15 


-2.91 + 0.03 


0.25 


32 


J082619. 70+314848 


08:26:19.70 


+31:48:48.0 


2.9122 


20.30+°!g 


< 12.60 


< 


-27.21 


35 


-1.88 + 0.15 


-2.01 + 0.01 


— 


32 


Q0836+11 


08:39:33.015 


+11:12:03.82 


2.4653 


20.581°;!° 


< 13.12 


< 


-26.97 


88 


-1.15 + 0.05 


-1.40 + 0.01 


0.59 


13 


J0929+2825 


09:29:14.49 


+28:25:29.1 


3.2627 


21.101°;°° 


13.15 ± 0.02 


-27 


47 + 0.02 


43 


-1.62 + 0.01 


-1.78 + 0.01 


— 


32 


BR0951-04 


09:53:55.69 


-05:04:18.5 


4.2029 


20.40+°;!° 


13.37 + 0.08 


-26. 


.54 + 0.13 


36 


-2.62 + 0.03 


< -2.57 


0.04 


6,13 


BRI0952-01 


09:55:00.10 


-01:30:06.94 


4.0244 


20.55_ 10 


13.55 ± 0.02 


-26 


52 + 0.10 


302 


-1.46 + 0.18 


-1.86 + 0.08 


— 


8,13 


PC0953+47 


09:56:25.2 


+47:34:44 


4.2442 


20.90+°-!g 


13.60 ± 0.10 


-26. 


.82 + 0.18 


70 


-2.19 + 0.03 


-2.52 + 0.08 


0.24 


22 


J1014+4300 


10:14:47.18 


+43:00:30.1 


2.9588 


20.50lHg 


12.76 ± 0.04 


-27 


25 + 0.04 


108 


-0.71 + 0.14 


-1.11 + 0.04 


— 


32 


Q1021+30 


10:21:56.84 


+30:01:31.3 


2.9489 


20.701°;!° 


< 12.91 


< 


-27.31 


100 


-1.94 + 0.02 


-2.16 + 0.01 


0.25 


13,22 


J103514.22+544040 


10:35:14.22 


+54:40:40.1 


2.6840 


20.501°;^° 


> 13.59 


> 


-26.42 


160 


-0.62 + 0.21 


-0.45 + 0.12 


1.43 


32 


Q1036-230 


10:39:09.4 


-23:13:26 


2.7775 


21.001°;!° 


12.97 + 0.05 


-27 


.55 + 0.11 


80 


-1.41 + 0.10 


-1.82 + 0.10 


— 


13 


Q1104-18 


11:06:32.96 


-18:21:09.82 


1.6614 


20.80+°,;!°, 


13.44 + 0.05 


-26 


88 + 0.11 


56 


-1.04 + 0.01 


-1.48 + 0.02 


— 


7 


BRI1108-07 


11:11:13.64 


-08:04:02.47 


3.6076 


20.501°;!° 


< 12.34 


< 


-27.67 


32 


-1.80 + 0.00 


-2.12 + 0.01 


0.19 


8,13 


J113130.41+604420 


11:31:30.41 


+60:44:20.7 


2.8760 


20.50+°'!* 


< 12.66 


< 


-27.36 


53 


-2.13 + 0.15 


-2.31 + 0.03 


— 


32 


HS1132+2243 


11:35:08.03 


+22:27:06.8 


2.7835 


2i.ooi°:°? 


< 12.69 


< 


-27.82 


56 


-1.99 + 0.07 


-2.22 + 0.02 


0.27 


22 


J115538. 60+053050 


11:55:38.60 


+05:30:50.6 


3.3268 


21.051°;!° 


13.73 ± 0.03 


-26 


84 + 0.10 


120 


-0.81 + 0.10 


-1.44 + 0.05 


1.21 


32 


Q1157+014 


11:59:44.81 


+01:12:07.1 


1.9440 


21.801°;!° 


> 14.80 


> 


-26.51 


84 


-1.36 + 0.06 


-1.81 + 0.04 


0.73 


9 


BR1202-07 


12:05:23.63 


-07:42:29.91 


4.3829 


20.601°;!^ 


< 13.06 


< 


-27.06 


170 


-1.81 + 0.02 


-2.19 + 0.12 


— 


2 


J120802. 65+630328 


12:08:02.65 


+63:03:28.7 


2.4439 


20.70+°'!* 


13.55 ± 0.03 


-26. 


.67 + 0.15 


60 


-2.32 + 0.15 


-2.55 + 0.01 


0.18 


32 


Q1215+33 


12:17:32.54 


+33:05:38.39 


1.9991 


20 95+° : ° 7 


< 13.17 


< 


-27.29 


42 


-1.48 + 0.03 


-1.70 + 0.05 


0.35 


6,13 


Q1223+17 


12:26:07.22 


+17:36:48.98 


2.4661 


21.501°;!° 


< 14.01 


< 


-27.01 


94 


-1.59 + 0.01 


-1.84 + 0.02 


0.62 


8,13 


Q1232+08 


12:34:37.55 


+07:58:40.5 


2.3371 


20.90+°,;!°, 


14.00 + 0.10 


-26. 


.42 + 0.14 


85 


-1.28 + 0.09 


-1.72 + 0.09 




9 


J1240+1455 


12:40:20.91 


+14:55:35.6 


3.1078 


21 30+"' 00 


> 14.34 


> 


-26.47 


335 


-0.85 + 0.03 


-2.11 + 0.02 




32 


J1240+1455 


12:40:20.91 


+14:55:35.6 


3.0241 


on 45+"- 00 


< 13.32 


< 


-26.65 


134 


-0.74 + 0.07 


-0.81 + 0.12 




32 


Q1331+17 


13:33:35.78 


+16:49:04.03 


1.7764 


21 14+ - 08 


< 13.54 


< 


-27.12 


72 


-1.42 + 0.00 


-2.02 + 0.00 


0.50 


6,13 



Quasar 


RA (2000) 


DEC (2000) 




Z ah., 


lo 


7#Ar- 


log 10 7V(CII*) C 


loglO^ 


Av 90 


[M /H] 




[Fc/H] 9 




W 


h 

L526 


Ref 


Q1337+11 


13:40:02.44 


+11:06:29.6 


2 


7959 


20 


95Io.!o 


13.11 ± 0.10 


-27 


36 ± 0.14 


60 


1.72 ± 


12 


—2.03 ± 


08 





25 


22,31 


BRI1346-03 


13:49:16.82 


-03:37:15.06 


3 


7358 


20 


_i_o.io 

'^-0.10 


12.55 ± 0.11 


-27 


69 ± 0.15 


38 


—2.33 ± 


01 


—2.63 ± 


02 





12 


6,13 


PKS1354-17 


13:57:06.07 


-17:44:01.9 


2 


7800 


20 


30-0.15 


12.76 ± 0.06 


-27 


06 ± 0.16 


30 


— 1.37 ± 


19 


— 1.79 ± 


05 





21 


22 


J141030+511113 


14:10:30.60 


+51:11:13.5 


2 


9642 


20 


°°-0.20 


< 13.01 


< 


—27.35 


46 


— 1.96 ± 


15 


—2.27 ± 


02 








J141030+511113 


14:10:30.60 


+51:11:13.5 


2 


9344 


20 


or, +0.15 
80 -0.15 


> 13.39 


> 


—26.92 


247 


—0.95 ± 


15 


—1.16 ± 


10 






32 


J141906. 32+592312 


14:19:06.32 


+59:23:12.3 


2 


2476 


20 


qc+0.20 
yo -0.20 


< 13.08 


< 


—27.39 


20 


—2.85 ± 
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are positive detections, green are upper limits, and blue 
are lower limits. The positive detections are distributed 



into two distinct regions centered near £ C =10~ 27A ergs 
s _1 H _1 and £ C =10~ 26S ergs s _1 H _1 . A histogram of 
the positive detections is plotted in Fig. O which shows 
the two peaks separated by a trough at £ c — ^ rlt RilO _27 ' 
ergs s _1 H _1 . Note, the trough is not an artifact caused 
by selection effects, since the signal-to-noise ratios of the 
C II* absorption profiles for the 5 DLAs with £ c ^£ c L rit 
are not exceptional for DLAs with positive detections. 
Standard statistical tests applied to these data support 
the bimodal hypothesis. We first used the KMM algo- 
rithm (Ashman et al. 1994; McLachlan & Basford 1987), 
which fits Gaussians to the data, computes the maximum 
likelihood estimates of their means and variances, and 
evaluates the improvement of that fit over the null hy- 
pothesis of a single Gaussian. We found the probability 
of the null hypothesis to be Pkmm(^c)=0.016; i.e., a sin- 
gle Gaussian fit can be excluded with 98.4 % confidence. 
This result is robust, as the value of Pkmm(4) is insen- 
sitive to a wide range of initial guesses for the means 
and variances of the input Gaussians. We then used the 
Bayesian mixture algorithm NMIX (Richardson & Green 
1997), which models the number of Gaussian components 
and mixture component parameters jointly and evaluates 
the statistical significance of these quantities based on 
their posterior probabilities. J. Strader (priv. comm.; 
2007) kindly analyzed the data with NMIX and found 
that the probability of the null hypothesis, Pnmix=0.08. 
The latter test yielded peak locations of 4=10~ 27 - 34±a06 
and 4=io- 26 - 58±0 06 ergs s" 1 H" 1 and values of 0.38 and 
0.62 for the sample fractions associated with the respec- 
tive peaks. 

These results were obtained by excluding both the 
lower limits and upper limits on £ c . This does not af- 
fect our conclusions provided the true values of £ c , i.e., 
£' rue , are drawn from the same parent population as the 
positive detections, which we assume is a bimodal distri- 
bution. Note by definition the values of ^* ruo equal the 
measured values of £ c in the case of positive detections, 
are less than the measured £ c for the upper limits, and 
are greater than the measured £ c for the lower limits. 
The distributions of lower limits and positive detections 
shown in Fig. [3^ are clearly compatible. However, be- 
cause the values of £* rue must exceed the corresponding 
lower limits, it is possible that the ^' ruo are drawn from 
a separate population concentrated above the peak at £ c 
w 10~ 266 ergs s _1 H _1 characterizing the positive detec- 
tions. But there is no sign in any of the other properties 
of DLAs with lower limits that distinguishes them from 
the DLAs with positive detections. This suggests that 
both sets of DLAs are drawn from a common popula- 
tion. In any case the small number of lower limits implies 
they have negligible impact on the case for bimodality or 
on the location of either peak inferred from the positive 
detections. 

By contrast, Fig. [3b shows that the £ r distribution of 



the upper limits peaks at 



10 



-27.0 



ergs s 



li- 



the 



same value of £ c where the positive detections exhibit a 
trough. The difference between the distributions is po- 
tentially important owing to the comparable sizes of the 
two samples. However, the distribution of upper limits 
shown in Fig. [3}d does not accurately represent the cor- 
responding distribution of ^' ruc . Because the values of 




-25 



log 10 1 [ergs s-' H"»] 



Fig. 3. — Histograms comparing positive detections (red) and 
(a) lower limits (blue), (b) upper limits (green), and (c) all the £ c 
values in Table 1 (magenta). 

£truo ghouid be l ower than the associated upper limits, 
the resulting distribution of £' rue should peak at values 
lower than shown in the figure. While we do not know 
how much lower, our recent (2006 Dec, 2007 April, and 
2007 Sept.) HIRES observations provide a clue. We ob- 
served 13 DLAs with upper limits set previously by data 
acquired with ESI. Using the higher spectral resolution 
of HIRES (i?=43,000 compared to i?=9000 for ESI), we 
detected four of these while the other nine remained un- 
detected. The HIRES observations moved five of eight 
ESI upper limits out of the trough to upper limits below 
the trough, while two of the four new detections were 
also moved below the trough and one remained within 
the trough. Therefore, it is plausible to assume that 
^true cou i(j typically be more than 0.4 dex lower than the 
upper limits; i.e., the evidence is consistent with a distri- 
bution resembling the positive detections by exhibiting a 
peak at 4«10 -27 - 4 ergs s" 1 H" 1 . 

Fig. [3J; plots the empirical £ c distribution for all 76 
DLAs in Table 1: the plot includes the positive detec- 
tions, lower limits, and upper limits. While the cor- 
responding £' ruo distribution for the entire sample has 
not fully been determined, the above arguments provide 
good evidence for a 'high cool' peak at £ c «10 _26 6 ergs 
H _1 and growing evidence for a 'low cool' peak at 
ergs s _1 H _1 ; i.e., a distribution resembling 



-27.4 



that of the positive detections. Although we do not know 
which values of £' rue to assign to the upper limits, we 
suggest that in most cases the range of values is given 



6 



by ^ rue =[10- 28 - 5 ,lCr 27 - ] ergs s' 1 H . Lower bounds 
on ^* ruc are set by the heating rates due to FUV and 
X-ray background radiation and excitation rates due to 
CMB radiation. The CMB sets a floor of ^ c ruc >10~ 28 
ergs s' 1 fT 1 for DLAs with z >3.5 (Wolfe et al. 2004; 
hereafter referred to as WHGPL), whereas X-ray heat- 
ing provides the lower bound for DLAs with lower red- 
shifts. We reject values of £ c lower than 10 -28 ' 5 ergs 
s' 1 H' 1 since the densities implied in the case of ther- 
mal balance n < 10~ 15 cm -3 (WHGPL). In that case 
the length scale d of gas clouds for DLAs with median 
column density Ahi=8x10 20 cm' 2 (WGP05), would ex- 
ceed 8 kpc. While the line-of-sight might traverse such 
distances through DLAs, the multi-component structure 
of the absorption-line profiles indicates the gas is con- 
fined to several smaller "clouds" characterized by a low 
volume filling factor (Nagamine ct al. 2007). This sce- 
nario is supported by multi-phase models in which the 
DLA gas density n > 10 cm' 3 and d < 30 pc (Wolfe, Ga- 
wiser, Prochaska 2003; hereafter referred to as WPG03). 
At the same time it is likely that £* ruc corresponding to 
most of the upper limits in Fig. [3J3 is less than 10 -27 ergs 
s _1 H _1 : since only four out of 32 upper limits exceeds 
2.Q-26.6 ergs s' 1 H _1 and we expect £' rue to be more than 
0.4 dex lower than the corresponding upper limit, then 
£true f or oyer qq% Q £ upper limits should be lower 
than 10" 27 ergs s" 1 H _1 . 

To summarize, standard statistical tests applied to the 
£ c distribution of positive detections support the bimodal 
hypothesis. The above arguments concerning the lower 
and upper limits are consistent with the hypothesis of 
two peaks in the £* ruc distribution. Because the pres- 
ence of two peaks is supported by physical arguments 
concerning the range of ^' rue , the case for bimodality is 
sufficiently compelling to consider independent tests to 
which we now turn. 

3. INDEPENDENT TESTS FOR BIMODALITY 

If the £ c distribution is bimodal, the critical cooling 
rate £ c * lt should divide the DLA sample into independent 



populations with physically distinct properties, 
section we investigate whether this is the case. 



In this 



3.1. Velocity Interval 

We first determine whether £" 11 divides DLAs into 
populations with distinct velocity structures. We choose 
the absorption-line velocity interval Augo as a measure of 
velocity structure, where A?;go is defined as the velocity 
interval containing the central 90% of the optical depth of 
unsaturated low-ion absorption lines (Prochaska & Wolfe 
1997): this definition guarantees that Ai>go corresponds 
to the bulk of the neutral gas producing the absorption. 
The measured values of Atigo are listed in column 8 of 
Table 1. To generate suitable subsamples we adopted 
the following criteria. For the positive detections we as- 
signed DLAs with £ c < £" lt to a 'low cool' subsample and 
DLAs with l c > £ c c Tit to a 'high cool' subsample. The di- 
vision is not as clear for DLAs with upper limits because 
£true ma y k e lower than l c J lt for objects with upper lim- 
its on £ c exceeding £" lt . To account for this possibility 
we assigned DLAs with upper limits given by logio^c < 
logio^c"' +Alogio^c rit to the 'low cool' subsample and 
those with higher upper limits to the 'high cool' subsam- 
ple. Similarly we assigned DLAs with lower limits given 
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Fig. 



- Histograms of velocity intervals Avgo divided such 
that the magenta histogram depicts DLAs with l c > l" 11 and the 
blue histogram depicts DLAs with t c < i^ Ilt . The inset shows the 
cumulative distributions for both subsamples. 

by logio4 < logio^c"' - Alog 10 ^ rlt to the 'low cool' sam- 
ple and those with higher lower limits to the 'high cool' 
subsample. The value of Alogio^™* is determined by the 
ratio of the £ c limit to our estimate of £' ruc . 
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Histograms with Alogi ^ rlt = 0.2 and ^ rit =10' 27 
ergs s' 1 H' 1 are shown in Fig. 01 In this case the 
'low cool' and 'high cool' histograms are clearly differ- 
ent. Whereas the bulk of the 'low cool' subsample clus- 
ters around At>g « 50 km s _1 , the bulk of the 'high 
cool' subsample DLAs exhibits a wider distribution cen- 
tered near 100 km s _1 . The difference between the two 
subsamples is also evident in the cumulative distribution 
shown in the inset in Fig. f5J The cumulative distribu- 
tions are obviously different with medians Awgo = 46 
km s' 1 for the 'low cool' subsample and 104 km s' 1 
for the 'high cool' subsample. Applying the standard 
Kolmogorov-Smirnov (KS) test we find the probability 
that the two distributions are drawn from the same par- 
ent population, PKs(Awgo) =2.7xl0 -5 ; i.e., the null hy- 
pothesis can be rejected at a high confidence level. 

To test the sensitivity of this conclusion to uncertain- 
ties in ££ rit and Alogio^"* we recomputed PKs(Avgo) 
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for a range of values compatible with the t c distribution 
in Fig. El i.e., 10~ 271 < £ c c Ht < 10~ 26 - 7 ergs s" 1 H" 1 
and < Alogi ^ rit < 0.4. The results shown in Table 2 
indicate that 8.5xlCT 6 < P K s(Av 90 ) < 1.5xl0~ 3 , with 
the largest value corresponding to the case Alogio^c" t= 0, 
which is unlikely. Next we checked whether the Avqo dis- 
tribution was bimodal with respect to other parameters 
measured in DLAs. The only parameter suitable for this 
purpose is redshift. Accordingly we divided our sample 
into low-z and high-z subsamples around the median red- 
shift z me( j=2.85. When we compared the resulting Avqq 
distributions, we found that Pks(Av9o) = 0.69. In other 
words velocity width is unlikely to be bimodal with re- 
spect to redshift. 
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Fig. 5. — Histograms of DLA metallicity divided such that the 
magenta histogram depicts DLAs with £ c > ££ rlt and the blue his- 
togram depicts DLAs with £ c < i^ xX . The inset shows the cumu- 
lative distributions for both subsamples. 



3.2. Metallicity 

Fig. [5] shows the metallicity distributions of the 'low 
cool' and 'high cool' subsamples. The logarithmic metal 
abundances with respect to solar, [M/H], are listed in col- 
umn 9 of Table 1, where M stands for an element found to 
be undepleted in DLAs such as S, Si, or Zn (Prochaska 
et al. 2003). As in Fig. |4] we divide the sample with 
^ rit =10- 27 ergs s- 1 R- 1 and Alogi ^ rit =0.2. Similar 
to the case of the Ai>go distributions we find the 'low cool' 
and 'high cool' distributions of metallicity are different. 
For these values of £^ rit and Alogio^"* the median metal- 
licity of the 'low cool' subsample is given by [M/H]i ou; = 
-1.74, and [M/H] w = -1.06 for the 'high cool' sub- 
sample. Applying the KS test, we find the probability 
that the two subsamples are drawn from the same parent 
population is given by Pks([M/H])=4.5x 10~ 5 . When 
we recomputed Pks([M/H]) for the range of ££ lt and 
Alogio^"* shown in Table 2, we found that 3.5xl0~ 6 < 
-Fks([M/H]) < 1.8xl0 -3 , where the largest value again 
corresponds to the unlikely case of Alogio^; nt =0. By 
contrast, when we split the sample according to red- 
shift we found that P K s([M/H])=0.50, indicating that 
DLA metallicities in the two redshift bins are not drawn 
from the same parent population. This may conflict 
with our detection of metallicity evolution with redshift 



(Prochaska et al. 2003) in which we found weak but sta- 
tistically significant redshift evolution; i.e., d[M/H]/dz — 
— 0.26±0.07. Part of the difference between the results 
could be due to sample size, since for the Prochaska et 
al. (2003) sample, which is about twice the size of the 
current sample, we find P K s([M/H])=0.001. While it 
is possible that some of the more recent measurements 
may dilute the original result, the vastly different values 
of Pks([M/H]) is puzzling. 
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Fig. 6. — Histograms of DLA dust-to-gas ratio divided such 
that the magenta histogram depicts DLAs with i c > and the 
blue histogram depicts DLAs with £ c < i c J l1 . The inset shows the 
cumulative distributions for both subsamples. 



3.3. Dust-to-Gas Ratio 

Not surprisingly the results for dust-to-gas ratio, K, re- 
semble the results for metallicity. This is shown in Fig.[6l 
which compares the k distributions for the 'high cool' 
and 'low cool' subsamples where k= lfj[ M / H ](10l Fo/Mlint - 
10 [Fc/M]) ( see WPG0 3) an d [Fe/M] and [Fe/M] int are 
the measured and intrinsic logarithmic abundances of Fe 
with respect to an undepleted a element, M: all abun- 
dances are with respect to solar ([Fe/H] is given in col- 
umn 10 in Table 1). We computed the numerical val- 
ues for k in Table 2 by assuming [Fe/M]; nt =—0.2 when 
[Fe/M] < -0.2 and [Fe/M] int =0 when -0.2 < [Fe/M] < 
0.0 (see WGP05). 

Fig. [6] shows that the dust-to-gas ratios of the 'high 
cool' subsample are significantly higher than for the 'low 
cool' subsample. We find Pks( k )= l.OxlO -3 for the dis- 
tributions shown in Fig. [BJ where £"' It =10 -27 ergs s _1 
tr 1 and Alogi ^ rit =0.2. We also find that 3.7xl0" 5 < 
Pks(k) < 1.4xl0~ 2 for the range in l c / lt and Alogi ^ rit 
in Table 2, with the largest values again given by the un- 
likely case Alogio^c rlt= 0. When we divided the sample 
according to redshift we found that Pks(«)=0.17. Al- 
though a unimodal distribution of k with respect to red- 
shift is less likely, it still cannot be ruled out with high 
significance. As a result bimodality is unlikely to arise 
with respect to redshift. 

3.4. Si II X 1526 Equivalent Width 

Prochaska et al. (2007) recently showed that DLAs 
exhibit a tight correlation between Si II A 1526 rest- 
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Histograms of Si II A1526 equivalent width divided 



such that the magenta histogram depicts DLAs with £ c > 



and 



the blue histogram depicts DLAs with £ c < £™ . The inset shows 
the cumulative distributions for both subsamples. 

frame equivalent width, W1526, and metallicity, [M/H]. 
Interestingly the correlation exhibits less scatter than 
the Avqq versus [M/H] correlation. This is a striking 
result because whereas Avqo and [M/H] are determined 
from the same neutral gas, W1526 is dominated by the 
kinematics of low optical-depth clouds that make an in- 
significant contribution to [M/H]. Prochaska et al. (2007) 
argue that W1526 is determined by the virialized random 
motions of clouds in the outer halo of the DLA galaxy, 
while Awgo is set by the motions (e.g. rotation) of the 
neutral ISM. As a result, the wide range of values ex- 
pected for impact parameter and galaxy inclination will 
cause sizable scatter in Avqq among galaxies with the 
same virial velocity. By contrast, the scatter in Wi 52e 
will be lower for sightlines encountering large numbers of 
randomly moving clouds in the halos of the same galax- 
ies. Because of the tight correlation between W1526 and 
[M/H], one might expect that W1526 determined for the 
'low cool' and 'high cool' DLAs are not drawn from the 
same parent population. 

Fig. [7] compares the W1526 distributions of the 'high 
cool' and 'low cool' subsamples for the standard val- 
ues ^rit =10 -27 ergs s -i H -i and Alog 10 ^ rlt =0.2. In 

this case the 51 measured values of W1526 are listed 
in column 11 of Table 1. The figure shows that 
the 'low cool' subsample exhibits systematically lower 
values of W1526 than the 'high cool' subsample: we 
find P K s(W^i526)=7\2xl0~ 4 . When we recomputed 
Pks(Wi 526 ) for the range of £f xt and Alogio^™* shown 
in Table 2, we found that 7.2xl0~ 4 < P K s{Wi 526 ) < 
6.5xl0~ 3 . When we split the sample according to red- 
shift we found Pks (W1526) =0.87, which further supports 
the hypothesis that the DLA sample is not bimodal with 
respect to redshift. 

3.5. H I Column Density 

Fig. [51 compares the Nm distributions of the 'high 
cool' and 'low cool' subsamples for the standard case 
4=10~ 27 ergs s" 1 H" 1 and Alogi ^ rit =0.2, where the 
values of Am are listed in column 5 of Table 1. By 
contrast with the previous results, the null hypothesis 
cannot be excluded at high confidence levels. Specifi- 
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Fig. 8. — Histograms of H I column density divided such that 
the magenta histogram depicts DLAs with £ c > £ c * lt and the blue 
histogram depicts DLAs with l c < £^. Tlt . The inset shows the cu- 
mulative distributions for both subsamples. 

cally, in this case Pks(Ahi)=0.87. More generally, 0.19 
< Pks(Ahi) < 0.96 for the range of values spanned by 
££ nt and Alogio^"* m Table 2. This has important im- 
plications to be discussed in § 5. When we split the data 
according to redshift we found that Pks(Ahi)=0.21, in- 
dicating that we cannot confidently rule out that Am 
values in the two redshift bins are drawn from the same 
parent population. As a result, the DLA sample is un- 
likely to be bimodal with respect to redshift. 

4. MODES OF HEATING 

In this section we describe the heating processes that 
balance cooling in DLAs. We interpret the physical sig- 
nificance of the critical cooling rate £" kt and discuss dif- 
ferences between the heat input into the 'low cool' and 
'high cool' DLAs. We consider two possibilities. 

4.1. Background Heating of 'Low Cool' DLAs 

The first possible explanation for the presence of a 
trough at £f h w 10" 27 ergs s" 1 H" 1 is that ^ rit is the 
maximum [C II] 158 (im emission rate of low-metallicity 
gas heated by background radiation at z ~ 3. Accord- 
ing to this scenario, DLAs with l c < ^ rlt , i.e., the 'low 
cool' DLAs, are neutral gas layers heated by background 
radiation alone, while DLAs with £ c > ^ rlt , the 'high 
cool' DLAs, are in addition heated by internal sources 
(see WHGPL and Wolfe & Chen 2006; hereafter referred 
to as WC06). 

This interpretation of ^ rlt is illustrated in Fig. [51 which 
shows £ c and the total heating rate, r to t, as functions 
of gas volume density n. We assume that F to t=rxR+ 
r pe +rc J where Txr, r po , and r c o are heating rates due 
to X-ray photoionization, grain photoelectric emission, 
and photoionization of C° (Wolfire et al. 1995). The 
curves in Fig. \§\ are solutions to the thermal balance and 
ionization equilibrium equations for neutral gas layers 
with Nm=10 20 - 6 cm~ 2 that are heated and ionized by 
background radiation fields computed for 2=2.3 (Haardt 
& Madau 1996: and more recently using CUBA 5 ). The 



5 C UBA (Ha ardt fc Madau 2003) is available 
http: / / pitto.mib.infn.it / ~haardt / cosmology.html 
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Fig. 9. — Thermal equilibria computed for neutral gas at z 
~ 3 heated by background radiation. Dominant heating mecha- 
nism is production of primary electrons by photoionization of H 
and He by soft X-rays {hv > 0.4 kev). Primary electron energy 
degraded by photoionization and collisional excitation processes, 
leading to the production of secondary electrons that heat ambient 
electrons by Coulomb interactions. Grain photoelectric heating by 
FUV radiation is also included. Background radiation fields com- 
puted by Haardt & Madau (2003; using CUBA). Panels show the 
background heating rates per H atom (solid black) and the [158] 
(im cooling rates, i c , per H atom (dotted blue) as functions of gas 
density. Panels show results for several values of gas metallicity, 
[M/H], Notice, £ c does not exceed 10~ 27 ergs s" 1 H _1 . See text 
for physical interpration of these curves. 



panels in this figure show the resulting r tot (n) and £ c (n) 
for metallicities spanning the range of values found in 
DLAs. 

Fig. [9] illustrates why background heating is an attrac- 
tive explanation for DLAs with £ c < £^ rlt . In particular £ c 
remains below 10 -27 ergs s _1 H _1 for the entire range of 
densities and metallicities in the figure. At densities be- 
low fXmaxj the lowest density where £ c is a local maximum 
(=£™ ax ), the gas is a warm neutral medium (WNM) in 
which T w 8000 K and £ c « r to t, since cooling is dom- 
inated by Lya and electron recombination onto grains. 
At densities above n max the gas is a cold neutral medium 
(CNM) in which T rs 100 K and £ c = T tot , because cooling 
is dominated by 158 /im emission. While a detailed dis- 
cussion about the behavior of £ c (n) is beyond the scope 
of this paper (see WHGPL and Wolfire et al. [1995] for 
more extensive discussions), we wish to emphasize two 
points. First, although £™ ax increases with increasing 
[M/H], Fig. [5] shows that it never exceeds 10 -27 ergs 
s _1 H _1 . Second, because £ c undergoes a sharp decline 



with decreasing density at n < n max , the range of 158 /im 
cooling rates available in the WNM for low values of £ c is 
much larger than the limited range of CNM cooling rates 
available at n > n max - Combined with the upper limits 
on density required in the case of background heating to 
explain the large observed ratios of C 11/ C I (WHGPL), 
this suggests that DLAs with £ c < ^ rlt are low-density 
WNM gas layers without local heat input. 

Because these conclusions are based on specific values 
of ./Vhi and z, we tested their generality by computing 
^max f or 40 DLAs in the standard 'low cool' popu- 
lation (see § 2). We adopted the metallicity, dust-to-gas 
ratio, and Ahi appropriate for each DLA, exposed the 
gas to the background radiation intensity inferred for the 
DLA redshift (Haardt k Madau 2003; using CUBA), and 
then computed £™ ax for the £ C M equilibrium curves. We 
found that £™ ax exceeded 10 -27 ergs s _1 H _1 for only 5 
of the 40 DLAs. This is consistent with our expectation, 
owing to the low metallicity of the 'low cool' population. 

However, despite the attractive features of this model 
it has serious problems. While it predicts £ c to be less 
than £™ ax , the measured values of £ c exceed €™ ax for 9 out 
of the 14 positive detections in the 'low cool' population. 
This indicates the presence of local heat sources for these 
DLAs, which contradicts the assumption of background 
heating alone. Of course this conclusion is based on our 
determinations of £™ ax which for low-metallicity gas is 
a sensitive function of the X-ray background intensity 
at z«3 and photon energies hv > 0.5 keV. Because the 
background intensity depends on volume emissivity, the 
critical quantity is the normalization of the X-ray lumi- 
nosity function for high- 2: AGNs at the characteristic lu- 
minosity, L*. Results from the recent study of Silverman 
et al. (2006; 2007) suggest that the Haardt-Madau (2003; 
using CUBA) determination, which was based on earlier 
work by Ueda et al. (2003), is a factor of 3 higher than 
the Silverman et al. (2006; 2007) determination for 3.0 
< z < 4.0, the relevant redshift interval for background 
radiation at z w 3. In addition the Haardt-Madau (2003; 
using CUBA) background must be increased by at least 
a factor of three for background heating to explain the 
measured values of £ c for the 'low cool' population. As a 
result, the X-ray luminosity function must be a factor of 
9 higher than the Silverman et al. (2006; 2007) value to 
explain the £ c values. While there are other uncertainties 
in determining the heating rates such as the carbon abun- 
dance and the grain-photoelectric heating efficiencies, we 
conclude that while it cannot be ruled out altogether, the 
background heating hypothesis is not likely to be correct. 



4.2. Local Heating 

We next consider the alternative hypothesis that the 
'low cool' DLAs are heated primarily by local sources. 
Suppose the background intensity, J„ kd , at z « 3 is half 
the Haardt & Madau (2003; using CUBA) intensity, in 
agreement with observational uncertainties. In that case 
the measured values of £ c exceed £™ ax for 11 of the 14 pos- 
itive detections, indicating that local energy input is re- 
quired for most of this population. To determine the level 
of input we adopted the model of WPG03 in which the 
star formation rate per unit area projected perpendicular 
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to a uniform gaseous disk, Esfr 6 , generates FUV radia- 
tion that heats the gas by the grain photoelectric mech- 
anism. Heat inputs by cosmic rays and locally generated 
X-rays are also included and are assumed to be propor- 
tional to Esfr- The § as ^ s a ^ w0 P nase medium in which 
the CNM and WNM are in pressure equilibrium at the 
pressure given by the geometric mean of the minimum 
and maximum pressures characterizing the equilibrium 
P(n) curve (Wolfire et al. 2003). Because the detected 
C II* absorption likely arises in the CNM, WPG03 were 
able to deduce unique values for the local FUV radiation 
intensity J^ 003,1 from measured values of £ c . Specifically, 
in the case of thermal balance WPG03 found that 



(2) 



where J v = J^ kd + J l ° c&[ and e is the photoelectric heat- 
ing efficiency (e.g. Bakes & Tielens 1994; Weingartner 
& Draine 2001). To obtain E SF r, WPG03 solved the 
transfer equation in a uniform disk and found that in 
the optically thin limit 
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where C (=8.4xl0- 16 ergs 

cm _2 s _1 Hz _1 [M Q yr _1 kpc~ 2 ]~ 1 ) is a conversion 
constant (see WC06), and ^dla and h are the radius 
and scale-height of the DLA disk. 

Applying the same technique to the DLA sample in 
Table 1 and assuming rDLA/^=50, we determined val- 
ues of Esfr for 38 of the 44 DLAs with positive detec- 
tions and lower limits. 7 The resulting distribution of 
Esfr is shown in Fig. [TUl where we have plotted the re- 
sults for the 'low cool' and 'high cool' populations sepa- 
rately. We include results for the 'high cool' DLAs, since 
background heating is insufficient to balance the 158 /im 
cooling rates from this population, indicating that local 
heating is also required in this case. The Esfr distribu- 
tion is of interest for several reasons. First, application 
of the KS test provides strong evidence for two distinct 
modes of star formation, since the probability that the 
two populations are drawn from the same parent distri- 
bution is given by Pks(Esfr)=3.4x10~ 7 . Furthermore, 
comparison of Fig. [10] with Figs. 3 to 6 shows evidence 
for less overlap between the 'low cool' and 'high cool' dis- 
tributions of Esfr than for Augo, [M/H], k, and IUi526- 
Indeed, comparison with Fig. 2 suggests that the Esfr 
and l c distributions are both bimodal. 

Second, although the WPG03 model assumes in situ 
star formation throughout the neutral-gas disk, the sur- 
vey of the Hubble Ultra Deep Field (UDF) by WC06 
places upper limits on such star formation, which rule 
out some of the scenarios discussed above. Specifi- 
cally, WC06 searched the UDF F606W image for low 

6 In WC06 we denoted this quantity by (ip*)±_ to distinguish it 
from i/>* , the SFR per unit area determined from the observed H I 
column density by assuming the Kennicutt-Schmidt relation. Since 
such distinctions are not relevant here and to be consistent with 
the more general usage (cf Kennicutt 1998) we adopt the Esfr 
notation. 

7 We could not determine Esfr for the remaining 6 DLAs either 
because background heating was sufficient to balance the observed 
cooling, or we could not deduce a dust-to-gas ratio from our deple- 
tion model since [Fe/H] > [M/H], 
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Fig. 10. — Histograms of EgpR obtained using the C II* model 
of WP03 and WGP03. Blue histogram refers to 'low cool' DLAs. 
Magenta histograms refer to 'high cool' DLAs, where solid and 
dotted histograms respectively depict lower limits and true values 
of of Esfr- 

surface-brightness emission from DLAs in the redshift 
interval z— [2.5,3.5] with angular (linear) diameters be- 
tween 6» DLA =0.25 " (d DLA =2 kpc) and 6» DLA =4.0 " 
(<^dla=31 kpc). They optimized the survey sensitivity 
by convolving the F606W image with Gaussian smooth- 
ing kernels with FWHM diameters #ker n = Odla'- the 
resulting «5cr surface brightness thresholds varied be- 
tween ^ rosh =28.0 mag arcsec- 2 for 6> kcrn =0.25" to 
/U thresh =2 9_7 mag ar csec- 2 for kem =4.O". They found 
the number of detected objects to decrease steeply with 
increasing #kem and to be at least two orders of magni- 
tude lower than predicted, which is determined by di- 
viding the large DLA area covering factor, C A =0.33 (for 
z=[2.5,3.5]), by the area per DLA, 7t#dla 2 /4. Stated 
differently, in situ star formation is unlikely to occur 
in DLAs at values of Esfr above the SFR thresholds, 
Egpjj , set by the UDF threshold surface brightnesses. 
To compute E^ sh we note that /„ =CEsfr[47t(1 + 
z) 3 cos(i)] _1 where I Vo is the intensity observed from a 
disk at redshift z, and i is the disk inclination angle. 
Therefore, a disk with a given Egp r ^ sh generates a range 
of intensities that depend on i. To assure the detection 
of a significant fraction of disks with threshold SFRs 
we set i=60°. This guarantees that I Vo > /^ rcsh (or 
HV < ^ rcsh ) for half the disks for which E SF r =E^ sh 
and i > 60°. 

The resulting Egp r ^ sh are plotted as functions of #kern m 
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Fig. 11. — Solid curves are threshold values of Esfr versus 
diameter of smoothing kernel used by the search of WC06. Results 
shown for z = 2.5 (blue), 3.0 (green), and 3.5 (red). Horizontal 
dotted lines depict cooling rates per H atom characterizing 'low 
cool' and 'high cool' populations. 

Fig. HU for three redshifts spanning the redshift interval 
of the WC06 survey. 8 The figure shows that S|^ sh de- 
creases in the interval #kern=[0.25",2.0"]. The decrease is 
caused by the Poisson error in surface brightness, which 
decreases with increasing aperture size. The decline of 
£gFjf h ceases at 0^ ein > 2.0" due to the appearance of 
the systematic error, which WC06 attribute to confusion 
noise arising from the low surface-brightness outer re- 
gions of bright galaxies. Comparison with Fig. [TOl shows 
that the SFR per unit area for the bulk of the 'high cool' 
population, exceeds 10~ 22 M Q yr _1 kpc~ 2 . Since this 
exceeds Egpjj 8 for all values of the redshifts and #kern 
depicted in Fig. 111! such star formation would have been 
detected by WC06. Because WC06 detected less than 
1% of the DLAs predicted for this population, we con- 
clude that the heat input from in situ star formation in 
'high cool' DLAs cannot balance the rate at which they 
cool. This is in agreement with the results of WC06 
who worked with the Kennicutt-Schmidt relation rather 
than the C II* technique. As a result, the FUV inten- 
sity that heats the gas is likely emitted by compact (i.e., 
#fwhm < 0.25 "), centrally located, regions that WC06 
identified as Lyman Break Galaxies (hereafter referred 
to as LBGs; Steidel et al. 2003). Note, in that case 

8 For a given surface brightness the above definition results in 
values of Egp r ^ sh that are a factor of 2.8 smaller than the effective 
SFR per unit area, tp*, discussed in WC06 (see their eq. 3). 



the Esfr obtained for this population (see Fig. [TOj) are 
lower limits on the true values of Ssfr, which we esti- 
mate from the following argument. Assuming the 'high 
cool' DLAs account for half the DLA population and 
contain LBGs, the ratio of areas occupied by DLAs to 
LBGs, 7r^DLA/ 7rr LBG' ^ s gi ven by the ratio of area cov- 
ering factors, Cpla COo1 /Clbg- Since WC06 found that 
Cdla/Clbg=330 in the redshift interval z= [2.5,3.5], we 

find that 7rr 2 )LA /7rr 2 jBG =165, where C^ ) 1 ^ :ool = 0.5Cdla- 
As a result, the true values of Ssfr are a factor of 165 
higher than the lower limits for the 'high cool' popula- 
tion. The distributions of both the lower limits and true 
values of Ssfr are depicted in Fig. [I0l 

Fig. [10] also shows that £ SFR > 1CT 3 M Q yr" 1 kpc -2 
for the bulk of the 'low cool' population. From Fig. [TT] 
we see that this exceeds £gp r R s h for 9^ eTn > 0.7" at z=2.5, 
0kem > 1.3" at z=3.0, and 6> kcrn > 2.0" at z=3.5. There- 
fore, in situ star formation is detectable from 'low cool' 
DLAs with angular diameters #dla above these limits. 
Since less than 1 % of the predicted number of DLAs 
were detected by WC06 and since the 'low cool' popula- 
tion comprises about half of all high z DLAs, we conclude 
that in situ star formation for this population is compat- 
ible with the observations only if the values of #dla are 
smaller than these limits. Note, we cannot rule out the 
possibility that similar to the 'high cool' population, star 
formation in 'low cool' DLAs is sequestered away from 
the DLA gas in compact LBG cores. However, owing 
to the different properties of the two populations, it is 
more likely that within the context of the local heating 
hypothesis the 'low cool' DLAs undergo in situ star for- 
mation. 
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TABLE 3 
Population Properties 


Property (median) 


Population 


'Low Cool' 'High Cool' 


logio4 [ergs s" 1 H~ x ] 
Av 90 [km s _1 ] 
[M/H] 
W 1526 [A] 
logioAfHl [cm -2 ] 


-27.29±0.07 
46±10 

-1.74±0.19 
0.26±0.09 
20.7±0.17 


-26.54±0.13 
104±15 

-1.06±0.13 
0.67±0.19 
20.7±0.19 


loglQK 

logi [4 ocal /10- 19 (cgs)] 
log 10 S S FR a [M Q yr- 1 kpc" 2 ] 


-2.57±0.17 
-18.58±0.17 
-2.78±0.34 


-1.60±0.26 
-17.71±0.12 
>-1.85 



5. SUMMARY OF DLA PROPERTIES 

At this point we pause to give an overview of the DLA 
properties discussed above. A summary is given in Table 
3 which lists the median values of the directly measur- 
able properties £ c , Av 90 , [M/H], Wi 52 e, and AT HI , and the 
model-dependent quantities k, J y ° cal , and Ssfr- With 
the exception of Ahi, all the properties of the 'high cool' 
DLAs have systematically larger values than their 'low 
cool' counterparts. Because Aiigo and W1526 measure 
spatial variations in velocity it is plausible to assume 
they measure mass. In that case the median halo mass 
of the 'high cool' DLAs would be 30 times larger than 
that of the 'low cool' DLAs. This follows from the as- 
sumptions that (1) Augo and W1526 are each proportional 
to the halo circular velocity w C i rc and (2) Mdm oc v^ ilc 
(see Mo et al. 1998) where Mdm is the halo dark-matter 
mass. A higher mass for the 'high cool' population is 
also consistent with its higher metallicity, since metallic- 
ity is an increasing function of mass in modern galaxies 
(Tremonti et al. 2004), and chemical evolution in ACDM 
models predict a similar trend at high redshifts (e.g. 
Nagamine 2002). Because the dust-to-gas ratio k traces 
metallicity, the higher median value of k for the 'high 
cool' DLAs is expected. On the other hand higher val- 
ues of £ c and Ssfr in the 'high cool' DLAs goes against 
the trend in current galaxies where star-formation activ- 
ity is more common in low-mass late type galaxies than 
in massive early- type galaxies (Kauffmann et al. 2003). 

While we have interpreted the two types of DLAs 
as physically distinct populations, let us consider the 
alternative hypothesis that they are distinguished by 
the difference in impact parameters of sightlines pass- 
ing through a unique population of DLAs. Suppose 'low 
cool' DLAs are metal-poor outer regions connected to 
metal-rich inner regions of the 'high cool' DLAs by a neg- 
ative metallicity gradient. In this case the lower values of 
Ai>go in the 'low cool' DLAs would arise from passage of 
QSO sightlines through neutral gas at larger impact pa- 
rameters where smaller changes in velocity generated by 
the systematic motions of the neutral gas are expected 
(Prochaska & Wolfe 1997). There are, however, at least 
two difficulties with this idea. First, since the values of 
W1526 are plausibly determined by the virialized random 
motions of clouds in the outer halo, W1526 would not be 
affected by a change in impact parameter, since tdla 
is small compared to the radius of the halo. Yet the 
median value of 1^1526 for the 'low cool' DLAs is is sig- 
nificantly lower than for the 'high cool' DLAs (see Table 
3). Secondly, one would also expect a decrease in Ahi to 



accompany the decrease in [M/H] with increasing impact 
parameter, which is contradicted by the good agreement 
between the median values of TVhi in Table 3 (but see 
Johansson & Efstathiou 2003). Third, it is difficult to 
understand how a difference in impact parameter could 
produce bimodality in £ c . 

Therefore, while we cannot rule out the alternative pos- 
sibility that the velocity fields in DLAs could be due to 
non-gravitational motions (e.g. galactic-scale winds, see 
discussion in § 7) , we shall adopt the working hypothesis 
that the crucial quantity distinguishing the two popula- 
tions is mass. 

6. CONNECTIONS TO BIMODALITY IN GALAXIES 

In this section we discuss analogies with bimodality in 
modern galaxies, and then place the results within the 
context of galaxy formation. 

6.1. Bivariate Distributions 

The bimodality of the £ c distribution in Fig. [2] brings 
to mind bimodal distributions in contemporary galax- 
ies. From their analysis of ^10 5 galaxies drawn from the 
SDSS survey, Baldry et al. (2004; see also Blanton et 
al. 2003) find the u — r color distribution to be bimodal. 
The distribution exhibits a blue peak at u — r rs 1.4 and 
a red peak at u — r ~ 2.4 . In a similar analysis, Kauff- 
mann et al. (2003) find bimodality in the distribution of 
the age parameter D n (4000) and the indicator of recent 
star-formation activity, H5a- The galaxies centered at 
the blue peak are roughly the same objects exhibiting 
peaks at young stellar ages and recent star-formation, 
while the the galaxies centered at the red peak exhibit 
peaks at older stellar ages and low levels of recent star- 
formation activity. 

A possible link to bimodality in DLAs follows from 
the bivariate galaxy distributions in the {u — r,M*) 9 , 
(L> n (4000),M*), and (R5 A , M*) planes. Specifically 
galaxies in the (u — r,M*) plane divide into separate 
blue and red sequences. Within each sequence, u — r 
is correlated with M* . Bimodality in the u — r distribu- 
tion arises from extracting the colors of galaxies within 
bins of constant M*. As M* increases from low values, 
where the color distribution is dominated by the blue 
peak, a red peak grows while the blue peak declines in 
strength until at the highest masses the distribution is 
dominated by the red peak. Similar behavior is seen in 
the (D„(4000),M») and (HSa, M») planes, except that 
US a is anti-correlated with M* . 

In Fig. [12] we plot DLA analogues to these bivariate 
galaxy distributions. Following the discussion in § 3.4 we 
substitute W1526 for M* because of evidence that W1526 
is an indicator of dark-matter mass, Mdm (Prochaska et 
al. 2007), and since M» should be proportional to Mdm- 
Fig. [T2"k shows an analogous and unambiguous correla- 
tion between [M/H] and logioW^i526- By analogy with the 
galaxy result the W1526 bin corresponding to low values 
of W 1526 (logioWi52 6 (A) < - 0.8) produces an [M/H] 
distribution which peaks at [M/H] ~ —2.3, while the bin 
corresponding to higher values, i.e., logioH / i526(A) > — 

9 Baldry et al. actually use the u — r,M r plane where M r is 
the Petrosian absolute red magnitude. But the tight correlation 
between M, and M r implies that M r is an accurate indicator of 
stellar mass. 
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0.2, gives rise to an [M/H] distribution, which peaks at 
[M/H] ~ —0.7. Although bins at intermediate values of 
W / i526 do not result in clear evidence for two peaks in the 
[M/H] distribution, this may be a consequence of small 
numbers statistics. Furthermore, small numbers statis- 
tics may be responsible for the absence of a clear discon- 
tinuity between the 'blue' and 'magenta' distributions in 
the ([M/H], logioWisze) plane. 
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Fig. 12. — Bivariatc distributions of various parameters versus 
Wis26- Plots are DLA analogues of bivariate distributions of var- 
ious galaxy parameters versus stellar mass, (a) Metallicity versus 
Wi526- Blue and magenta points are 'low cool' and 'high cool' 
populations with measured values of W1526 ■ Note the strong posi- 
tive correlation between [M/H] and W15261 and the dominance of 
'low cool' points for low values of W1526 and 'high cool' points 
for high values of Wi526- (b) [f c /ergs s — 1 H — I ] versus Wi526- In 
this case only DLAs with positive detections of C II* are shown. 
The discontinuity between the two populations is obvious but no 
evidence for correlations is present, (c) Psfr/M© yr — 1 kpc — 2 ] 
versus VFl526- Systems with high Ssfr dominate at large values 
of Wi$26- Note, values of Esfr for 'high cool' DLAs are lower 
limits. 

On the other hand such a discontinuity is present in 
Fig. [T2"b. which plots logio^c versus logioWi526- In this 
case we only plot positive detections, since the upper 
and lower limits would obscure possible trends in the 
(£ c ,Wi526) plane. While neither population shows evi- 
dence for correlations between these variables, the bin 
with logioW / i526(A) > —0.2 produces an l c disribution 
dominated by large cooling rates. Interestingly, the 
absence of a clear anti-correlation between logio^c and 
logi H^i526 runs counter to the trend in the (HSa, M*) 
plane, which is the closest galaxy analogue, as we have 
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Fig. 13. — Same as Fig. 1121 except that Augo is substituted for 
W1526 as the independent variable. 

interpreted £ c as an indicator of star formation. Fig. 112b 
plots logioSsFR versus logio 1^1526, which is also analo- 
gous to the bivariate distribution in the (H^, M*) plane 
since Ssfr is, of course, a star-formation indicator. In 
this case, evidence for discontinuities between the two 
populations is strengthened by arguments in § 4.2, which 
imply that the values of Esfr for the 'high cool' DLAs 
in Fig. [12] are conservative lower limits. We find similar 
results, with somewhat larger scatter, when Avgo is used 
as a mass proxy, as shown in Fig. 1131 

However, we wish to point out an important difference 
between the C II* and traditional star-formation signa- 
tures. While the HSa index is sensitive to star bursts 
ending ~ 10 s to 10 9 years before the detection of US ab- 
sorption in galaxies, the £ c cooling rate is sensitive only 
to star formation contemporaneous with the epoch at 
which C II* absorption is detected. The reason for this 
stems from the short cooling times in DLAs, where 



^cool — 



5kT 



(4) 



in the case of CNM gas. Here we have assumed the gas 
to be in pressure equilibrium (Wolfirc et al. 1995) and we 
have ignored the contribution of CMB excitations to the 
158 //m emission rate defined in Eq. 1 (WGP03). Assum- 
ing T = 100K, we find cooling times of « 3xl0 6 yr and 
2xl0 5 yr for the median 'low cool' and 'high cool' DLAs 
in Table 3. Therefore, our observations are unlikely to 
detect 158 /mi cooling rates from star bursts with dura- 
tion At short compared to the ~ 2xl0 9 yr time interval 
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corresponding to the rcdshift search window used to find 
the DLAs in Table 1. The implication is that while the 
star formation history of DLAs may be punctuated with 
isolated, short-lived starbursts, the C II* technique is 
sensitive only to a continuous mode of star formation, or 
multiple bursts with a high duty cycle. Consequently, we 
may have detected the latter modes of star formation. 

The short cooling times raise the possibility that the 
'low cool' and 'high cool' DLAs represent different lev- 
els of star formation activity in the same object. In this 
scenario an underlying base level of activity gives rise 
to the 'low cool' DLAs, while bursts in star-formation 
rates produce the 'high cool' DLAs. The l c values would 
track the star formation rates because of the short cool- 
ing times. However, the metallicity distributions of the 
two populations would be indistinguishable since metal- 
licity is a byproduct of star formation history rather than 
instantaneous star formation rate. While the star bursts 
in the 'high cool' mode could affect the velocity structure 
of the gas through stellar winds, the high duty cycle of 
the bursts implies that the velocity structure of the gas 
is also a function of star-formation history rather than 
instantaneous rate. For these reasons it is difficult to 
understand how varying star formation rates in a sin- 
gle class of DLAs could explain why the metallicity and 
velocity structure of the two populations are so different. 

Therefore, in common with modern galaxies the bi- 
modality in DLAs may arise from the presence of two 
distinct sequences in mass and star-formation activity. 
The obvious question is whether the DLA phenomenon 
is the precursor to the bimodality detected in galaxies? 
Stated differently, have we detected the early stages of 
galaxy bimodality in DLAs? 

6.2. Galaxy Formation Models 

To answer this question we turn to models for the ori- 
gin of galaxy bimodality and then suggest an empirical 
test. 

6.2.1. Transition Mass Models 

Dckcl & Birnboim (2005) interpret galaxy bimodality 
in terms of a transition occurring at the 'shock' dark- 
matter mass, M s hock~10 115 M Q (Keres et al. 2005; 
Dekel & Birnboim 2005), corresponding to M* «10 10 5 
M Q , which is the stellar mass separating the two se- 
quences in the (D„(4000),M») and (H5a, M*) planes 
(Kauffmann et al. 2003). At high redshifts, halos with 
Mdm < Mshock accrete gas in a cold flow in which the 
temperature never climbs to the virial temperature of 
the halo: such processes result in the formation of disks 
observed along the lower-mass blue sequence. By con- 
trast, gas accreting onto high-redshift halos with Mdm 
> M s hock is first heated to the virial temperature (~10 6 
K) by stable shocks. The gas later contracts in a quasi- 
static spherically symmetric cooling flow to form the 
bulges found in the more massive galaxies populating 
the red sequence. In this case outlying cold gas accretes 
along filaments connecting the massive halos to the IGM. 

This scenario provides plausible explanations for sev- 
eral phenomena related to bimodality in DLAs. First, 
the cold gas accreting onto halos with Mdm < M s h oc k 
naturally evolves into disk-like structures, which act as 
neutral-gas reservoirs for in situ star formation. This is 
in accord with our finding that in situ star formation in 



DLAs is the dominant mode in the 'low cool' population, 
which we identify with halos having Mdm < M s h oc k- The 
low level of star formation predicted for these objects is 
probably related to the low molecular content of the DLA 
gas or to the increase with redshift of the Toomre critical 
surface density (see WC06). Second, the transition-mass 
scenario predicts little, if any, in situ star formation in 
the hot gas accreting onto halos with Mdm > M s h oc k- 
Rather, star formation in these objects is predicted to 
occur in dense, centrally located bulges that are fueled 
by filamentary cold flows penetrating the hot gas (Dckcl 
& Birnboim 2005; Keres et al. 2005). This fits in nat- 
urally with our finding that star formation in the 'high 
cool' population, which we identify with halos having 
Mdm > M s hock, mainly occurs in compact star- forming 
bulges sequestered away from the DLA gas detected in 
absorption, which we associate with the dense filamen- 
tary gas predicted in the simulations. This crucial result 
suggests that the red-blue bimodality in galaxies and the 
'high cool'/'low cool' bimodality in DLAs have a common 
physical origin. As a result bimodality in galaxies may 
originate in DLAs since it predicts that at large redshifts 
the most active star-forming objects are bulges located in 
the most massive halos. This is in contrast to bimodality 
in modern galaxies, in which the bulges and spheroids of 
massive galaxies exhibit little evidence of star formation. 
Presumably feedback by AGNs and supernovae that form 
in the bulges inhibits star formation subsequent to the 
epochs during which the outlying neutral gas is detected 
as high-redshift DLAs (e.g. Governato et al. 2004). 

If bimodality in DLAs stems from a transition in star- 
formation modes at Mdm~ M s hockj then about half of all 
DLAs would need to arise in halos with Mdm > M s h oc k: 
this conclusion follows from our finding that about half 
of our DLA sample is in the 'high cool' population. Be- 
cause of the large value of M s h oc k, the high median mass 
M me d predicted for DLA halos would seem to contra- 
dict the standard hierarchical paradigm in which M moc i « 
10 10 Mq (Haehnelt et al. 1998; Johansson & Efstathiou 
2006). However, the results of recent numerical simu- 
lations reveal a more complex picture. Nagamine et al. 
(2004) find that M mcc i is a sensitive function of feedback. 
The model in best accord with observations is their Q5 
run in which high- velocity (~ 400 km s _1 ) winds effec- 
tively eject most of the neutral gas from low-mass halos. 
In this case M mod « 10 115 M (Nagamine et al. 2007), 
which agrees with the values predicted for M s h oc k and 
is consistent with the range of masses deduced from the 
cross-correlation between DLAs and LBGs (Cooke et al. 
2006). While this model does not explain all the DLA 
data, in particular the predicted area covering factor for z 
= [2.5,3.5] is about a factor of 2 lower than observed, the 
overall agreement is sufficient to imply that such large 
values of M me d are plausible. 

6.2.2. O VI Absorption Tests 

An important feature of the transition-mass model is 
the presence of hot gas behind the accretion shock in ha- 
los with Mdm > M s h oc k- Hot gas is predicted only in 
this mass range, as halos with Mdm < M s h oc k accrete 
gas in cold flows even when heating due to supernova 
feedback is included (Keres et al. 2005). Consequently, 
hot gas may be a signature of the 'high cool' popula- 
tion. The recent detection of O VI absorption in a large 
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Fig. 14. — Comparison between [M/H] histograms for DLAs in 
'high cool' sample (green) with blue histogram depicting 9 DLAs 
reported to exhibit intervening O VI absorption and red striped 
subsample showing the 6 of these with robust detections of O VI 
absorption. 

fraction of DLAs by Fox et al. (2007) and their inter- 
pretation that O VI is collisionally ionized suggests that 
hot gas is present in many DLAs. Although Fox et al. 
(2007) suggest this gas is shock heated by supernova rem- 
nants, the upper limits on star formation in DLAs set by 
WC06 restrict this form of energy input, provided the 
WC06 limits apply to gas surrounding compact LBGs. 
While such limits do not yet exist, let us assume that O 
VI absorption in DLAs arises in gas heated by accretion 
shocks. In that case the DLAs with O VI absorption 
would belong to the 'high-cool' population. Searches for 
C II* absorption have been carried out for three of the 
six DLAs in which O VI absorption was unambiguously 
detected in intervening systems. In DLA0112— 306 at 
z=2.702 (Srianand et al. 2005), £ c exceeds £f it , while in 
DLA2138-444 at z = 2.852 (Srianand et al. 2005) and 
DLA2206-199 at z = 2.076 (see Table 1) £ c is less than 
£" lt . Thus, the C II* data do not provide a statistically 
meaningful test. 

We have considered an alternative test in which we 
compared the [M/H] distributions of DLAs in the 'high 
cool' population with those exhibiting O VI absorption. 
The results are shown in Fig.Qj] When all nine interven- 
ing systems are included in the O VI sample, the KS tests 
yields P K s([M/H])=0.009 for the probability that the two 
distributions are drawn from the same parent popula- 
tion. This suggests that O VI absorption need not occur 
in 'high cool' DLAs. On the other hand, three of the O 
VI identifications are problematic: in each case only one 
member of the O IVAA 1031, 1033 doublet was identified, 
possibly due to blending of the missing O VI line with 
Lya forest absorption lines. But, it is equally plausible 
that all three features were misidentified as O VI. In that 
case only six DLAs would exhibit O VI absorption (see 
Fig. EI}, and a revised KS tests yields P KS ([M/H])=0.15. 
As a result, the data are not inconsistent with the hy- 
pothesis that the 'high cool' DLAs are in massive halos. 

7. CONCLUSIONS 

In this paper we searched for evidence of bimodality in 
DLAs. Our principal diagnostic tool was the [C II] 158 



/im cooling rate per H atom, £ c , which we obtained from 
measurements of the C II* A 1335.7 and damped Lya 
absorption lines. In addition to the C II* transition we 
used accurate velocity profiles of resonance transitions to 
measure velocity widths, equivalent widths, abundances, 
and dust-to-gas ratios for 76 DLAs. Our results are sum- 
marized as follows. 

(1) Our studies of the C II* A 1335.7 absorption line 
resulted in 37 positive detections, 7 lower limits, and 32 
upper limits of £ c . The positive detections show strong 
evidence for a bimodal distribution with two peaks at 
4=10~ 27 4 and 10" 26 - 6 ergs s" 1 H" 1 separated by a 
trough at ^ rit «10" 27 ergs s" 1 H _1 . In § 2 we ar- 
gued that the distribution of the true values of l c corre- 
sponding to the lower and upper limits is consistent with 
the distribution of positive detections. Stated differently, 
there is compelling evidence that the £ c distribution in 
Fig. [2] is a faithful representation of the parent popula- 
tion from which all the true values of £ c are selected. 

(2) In § 3 we tested the bimodality hypothesis by the 
following method. We first divided the full DLA sam- 
ple into a 'low cool' subsample with £ c < £" lt and a 
'high cool' subsample with £ c > £" lt . We then compared 
distributions of various physical parameters in the two 
subsamples. The probability that the parameters Ai>go, 
[M/H], k, and W1526 in the two subsamples are drawn 
from the same parent populations is small, ranging be- 
tween 10~ 5 to «5xl0~ 3 . On the other hand the dis- 
tributions of A^hi are consistent with being drawn from 
the same parent population. Therefore, the two subsam- 
ples are likely to be separate populations with distinct 
distributions of velocity width, metal abundances, and 
gas-to-dust ratio, but similar H I column-density distri- 
butions. 

(3) In § 4 we considered different physical processes re- 
sponsible for heating the 'low cool' and 'high cool' DLAs. 
We studied the possibility that the 'low cool' DLAs were 
comprised of warm neutral-medium gas heated by back- 
ground radiation alone, but considered it to be unlikely 
because background heating does not provide sufficient 
energy input to account for the cooling rates in the ma- 
jority of these DLAs. 

(4) We then investigated whether local heating could 
account for the cooling rates in both DLA populations. 
Applying the two-phase model of WPG03, in which FUV 
radiation emitted by massive stars heats the gas by the 
grain photoelectric mechanism and assuming that C II* 
absorption arises in CNM gas, we found the resulting 
SFRs per unit area, Esfr, to be bimodal. We further 
found that the infrequent occurrence of extended low 
surface-brightness galaxies in the UDF (WC06) rules out 
in situ star formation as the heat source for the 'high 
cool' DLAs. Rather, these DLAs are likely heated by 
FUV radiation emitted by compact bulge sources, iden- 
tified as Lyman Break Galaxies, embedded in the neutral 
gas. On the other hand in situ star formation in the 'low 
cool' DLAs is compatible with the UDF results, provided 
the DLA diameters are less than ~ 10 kpc. The short 
cooling times of the gas further imply that the star for- 
mation we detect occurs continuously rather than in a 
few isolated bursts. We conclude that star formation oc- 
curs in all DLAs and proceeds in one of two modes, in 
situ or bulge dominated. 
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(5) In § 5 we compared the median values of proper- 
ties in the two populations. The results suggest that the 
crucial parameter distinguishing the two populations is 
mass. That is, 'high cool' DLAs are embedded in ha- 
los that are significantly more massive than the halos 
encompassing the 'low cool' DLAs. 

(6) The idea of mass sequences brings up analogies 
with bivariate distributions in modern galaxies. Recent 
surveys reveal separate parallel 'blue' and 'red' galaxy 
sequences in which galaxy color and age are correlated 
with stellar mass, while star-formation activity is anti- 
correlated with stellar mass. When the galaxy sample is 
divided into bins of stellar mass, the color distribution in 
each mass bin is bimodal, with an increase in the ratio 
of 'red' to 'blue' galaxies with stellar mass. In § 6 we 
constructed bivariate distributions for analogous proper- 
ties of DLAs. Using W1526 (or Avgo) as a proxy for mass 
and [M/H], as a proxy for u — r color, we found strong 
evidence for a correlation between [M/H] and Wi526- We 
also found tentative evidence for two parallel sequences 
in the (^ C) W^i526) an d (Esfr,Wi526) planes, where l c and 
Esfr are possible proxies for IWa, an indicator of recent 
star formation in galaxies. We found tentative evidence 
for (1) bimodal distributions of i c and SgFR in bins of 
fixed W1526, and (2) a larger fraction of 'high cool' DLAs 
with high values of £ c and Esfr in bins with the largest 
values of Wi526- Therefore, star-formation activity in 
DLAs may increase with mass, which is consistent with 
periods of high-z star formation in early-type galaxies. 

(7) In § 6 we placed these results in the context of 
current galaxy formation theory. We found that the 
transition-mass scenario, introduced to explain bimodal- 
ity in modern galaxies, provides a plausible scenario for 
the onset of bimodality in high-z DLAs. In this picture 
the 'high cool' DLAs comprise outlying neutral gas fila- 
ments penetrating hot, virializcd gas that fills halos with 
mass, Mum > M shock (=10 n - 5 M Q ). The inflow of these 
cold streams through the hot gas results in active star 
formation in compact bulges, which we presume are the 
embedded LBGs that heat the surrounding neutral gas. 
On the other hand accretion of gas onto halos with Mdm 
< M s hock only produces cold inflows that result in the 
formation of neutral disks. We associate these objects 
with the 'low cool' DLAs in which the gas is heated by 
in situ star formation. Since independent simulations 
demonstrate that half of the DLA population could have 
Mdm > Mhock, this scenario is in accord with the l c 
statistics. 

To summarize, the most significant result of this paper 
is that the bimodality observed in modern galaxies may 
originate in high-z DLAs. According to Dekel & Birn- 
boim (2006) at high redshift the modes of star formation 
bifurcate around halos with Mdm ~ M s h oc k: halos with 
Mdm > M s hock accrete gas in hot spherical inflows that 
are penetrated by cold streams generating starbursts in 
the inner disks or bulges, while the inflow of cold gas 
onto halos with Mdm < M s h oc k leads to disk growth 
and in situ star formation throughout the disks. This 
results in a bifurcation in cooling rates, £ c , for the fol- 
lowing reasons. First, the bifurcation of star formation 
modes around M s h oc k creates a bifurcation in < J l ° cai >, 
the locally generated mean intensity averaged over all im- 
pact parameters, which is given by the same expression 
for bulge dominated and in situ star formation modes 



(WGP03). In both cases < J l ° cal > oc M,/7rr DLA (Eq. 3), 
where M* is the total SFR and tola is the radius of the 
neutral-gas disk. While both M* and tola presumably 
depend on Mdm (Springel & Hernquist 2003; Mo et al. 
1998), they also depend on environmental factors, such as 
gas density, which clearly differ for the two populations. 
As a result, < J l ° cal > will not be a continuous function 
of Mdm near M s h oc k, but rather will bifurcate. Because 
£ c oc neJ l ° caX for photoelectrically heated gas in thermal 
balance (eq. 2) and the product ne is about the same in 
the 'high cool' and 'low cool' DLAs, the £ c distribution 
should exhibit an analogous bifurcation or bimodality. 
The transition to bimodality in modern galaxies occurs 
at lower redshifts, say z < 2, where feedback processes 
quench star formation at Mdm > M s hock (Dckcl & Birn- 
boim 2006) and the 'high cool' DLAs shut down their 
cooling rates and evolve onto the red sequence. 

Of course, there are caveats to these conclusions. First, 
the analogies with mass correlations in modern galaxies 
are tentative due to the limited number of DLAs with 
positive detections of C II* absorption. We need to in- 
crease the size of the sample to establish the reality of 
this result. 

The second caveat is related to our assumption that 
Awgo and VF1526 are measures of virial velocities, hence 
of dark-matter mass. Razoumov et al. (2007) describe 
the difficulties in constructing self-consistent numerical 
ACDM models in which the DLA velocity fields are 
caused by virial motions alone: the challenge is to re- 
strict DLA gas to halos sufficiently massive to generate 
significant virial velocities and sufficiently numerous to 
reproduce the observed DLA covering factors. An al- 
ternative explanation is that gas motions in DLAs are 
due to outflows; i.e. winds (see Nulsen et al. 1998; 
Schaye 2001). The presence of winds is plausible in the 
case of 'high cool' DLAs, since we have associated these 
DLAs with LBGs, which exhibit P Cygni profiles (Pet- 
tini et al. 2002). It is reasonable to assume that DLAs 
with higher SFRs would produce higher outflow veloci- 
ties, which would be consistent with the larger values of 
Avg and W1526 in 'high cool' DLAs. This possibility 
is also consistent with observations of nearby starburst 
galaxies, which exhibit correlations between outflow ve- 
locities of neutral gas and SFRs (Martin 2005). As a 
result, the parameters Awg and IU1526 could be direct 
signatures of outflows rather than virial velocities. How- 
ever, since starburst galaxies also exhibit a correlation 
between outflow velocity and circular velocity (Martin 
2005), the DLA parameters could be indirect signatures 
of dark-matter mass even in this scenario. Therefore, 
while the physical origin of the velocities within DLAs 
remains controversial, a bimodal distribution of dark- 
matter mass is the most plausible explanation for the 
bimodal properties of DLAs. 

We finish with the following question: where are the 
DLAs in which C II* absorption arises in WNM gas? In 
§ 4.2 we argued that C II* absorption in 'high cool' and 
'low cool' DLAs arises in CNM gas exposed, respectively, 
to high and low values of J y ° cal . But if DLAs are multi- 
phase media in which the CNM and WNM are in pressure 
equilibrium, we would expect a similar fraction of DLA 
sightlines to intersect each phase if, by analogy with the 
Galaxy ISM, they have comparable area covering factors 



17 



(McKee & Ostriker 1977). Because the grain photoelec- 
tric heating efficiency, e (Eq. 2), is considerably lower 
in the WNM than in the CNM (Weingartner & Draine 
2001), it is also possible for C II* absorption resulting in 
£ c < t c * lt to form in sightlincs intersecting only WNM gas 
in 'high cool' DLAs. The problem with this scenario is 
that it would also predict velocity widths, metallicities, 
and dust-to-gas ratios significantly larger than observed 
for the 'low cool' DLAs. Thus, while some DLAs with 
£ c < may arise in the 'high cool' population, most 
of them plausibly form in a separate population of CNM 
gas. Therefore, if the area covering factors of the WNM 
gas in both populations are significant, the heating effi- 
ciencies must be sufficiently low for the resulting values 
of £ c to be less than the 10 -27 4 ergs s _1 H _1 peak of the 
'low cool' population. Consequently, the answer to our 
question may be that the sightlincs encountering only 



WNM gas may be responsible for many of the upper 
limits in Fig. [TJ 
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